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INTRODUCTION 


This  report  is  a  summary  of  the  principal  research  findings 
completed  under  grants  N00014-95-1-0253  (February  2  1995  to  January 
31,  1996)  and  its  follow-on,  N000 14-96- 1-1 136  (August  1  1996  to  June  30, 
1997),  on  experimental  and  theoretical  research  on  Free  Electron  Lasers. 

In  this  two  and  one-half  year  period,  work  was  completed  on  a  thesis 
project  involving  the  generation  of  FEL  harmonic  radiation,  and  additional 
work  was  done  on  the  theory  of  short-pulse  operation  of  the  FEL.  A  more 
recent  study,  involving  the  production  of  Cerenkov  radiation  from  short 
pulses  of  high  energy  electrons,  was  partly  supported  by  this  research  as 
well.  The  major  accomplishments  of  this  projects  are  given  here  as 
reprints  in  the  next  three  sections. 

Section  I  consists  of  a  paper  showing  how  a  50psec  pulse  of  high 
intensity  millimeter  radiation  (wavelength  =  1.5mm,  3GW/cm2)  could  be 
produced  from  a  750kV  electron  beam;  it  was  presented  in  SPIE  section  on 
intense  microwave  pulses  at  San  Diego  in  July  1995,  and  is  published  in  the 
SPIE  Volume  2555,  p.  250.  This  work  is  an  outgrowth  of  work  done  on 
short-pulse  FEL  amplifiers  that  was  begun  with  support  from  the 
preceding  ONR  grant,  and  which  is  described  in  publications  reproduced  in 
the  final  technical  report  (Columbia  Plasma  Laboratory  report  #128, 

1995). 


Section  II  is  a  preprint  of  a  paper  that  has  been  accepted  for 
publication  in  the  IEEE  Transactions  on  Plasma  Science,  special  issue  on 
high  power  microwave  generation  (to  be  published,  1998).  The  paper 
describes  a  source  which  can  generate  MW  level  pulses  of  psec  broadband 
microwave  radiation,  using  a  pulsed  electron  beam  having  6MeV  energy. 

Section  III  is  a  reprint  of  a  paper  published  in  Physical  Review  E 
summarizing  our  findings  on  a  special  type  of  harmonic  FEL  radiation;  it 
also  appeared  in  earlier  versions  in  the  Nuclear  Instruments  and  Methods 
volumes  relating  to  the  FEL  meetings  in  the  summers  of  1996  and  1995 
(Rome,  and  New  York  City  meetings  respectively),  as  well  as  in  the  Beams 
‘96  Conference  held  in  Prague.  We  reproduce  the  Physical  Review  paper, 
as  it  is  the  latest  and  most  comprehensive  version.  This  work  is  an 
outgrowth  of  our  interest  in  FEL  harmonic  radiations  which  dates  from 
1991.  We  use  a  “waveguide  FEL”  which  is  driven  by  a  magnetron  at  its 
lower  frequency  “intersection”  (unstable  root)  at  24Ghz,  and  which  is 
configured  as  a  travelling  wave  amplifier.  This  radiation  bunches  the  beam 
at  this  frequency,  and  its  harmonics.  The  bunch  harmonics  are  arranged  to 
couple  to  the  third  harmonic  resonant  FEL  interaction  (high  frequency 


“intersection”)  at  72Ghz,  and  we  find  a  certain  amount  of  power  there  that 
is  coherent  with  the  low  frequency  source.  A  more  unusual  finding  is  the 
non-resonant  excitation  of  the  seventh  harmonic  of  the  24Ghz  source, 
which  generated  appreciable  amounts  of  power.  This  occurs  because  the 
bunched  beam  harmonic  source  currents  travel  at  the  same  speed  as  the 

seventh  harmonic  TE72  wave . an  artifact  of  cylindrical  geometry  and 

Bessel  function  roots.  To  excite  the  seventh  harmonic,  the  electron  beam 
must  be  displaced  off-axis.  These  phenomena  were  studied  using  a  ID  FEL 
slippage  two-frequency  theory;  the  results  were  checked  by  others  using  a 
3D  theory.  Thus  we  have  demonstrated  new  mechanisms  for  generating 
coherent,  phase-stable  radiations  at  millimeter  wavelengths  using  the  FEL. 
It  is  worth  noting,  that  another  group  has  recently  reported  very  tight 
bunching  of  a  beam  at  the  lower  intersection  frequency,  with  substantial 
harmonic  content  [Physical  Review  Letters  79,  3905  (1997)]. 

* 

The  project  supported  the  part-time  work  of  a  visiting  scholar  (Dr. 
T-B.  Zhang),  Professor  T.C.  Marshall,  and  a  full-time  graduate  research 
assistant,  who  was  graduated  with  the  PhD.  in  June,  1997. 
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Narrow  high  power  microwave  pulses  from  a  free  electron  laser 

T.  C.  Marshall  and  T.  B.  Zhang 

Department  of  Applied  Physics,  Columbia  University 
New  York,  New  York  10027 


ABSTRACT 

We  have  explored  high  power  microwave  (X  =  1  5mm ")  nul v  i 

the  ID  Compton  EEL  equations  with  slippage.  For  an  appropriate  moer  JarfSSv  ^T-r*  apered  undulator  PEL  using 

°a"?*an  Propagaie  in  a  nearly  ^iSniTwa,  £Vt  ^s“  s£™,^  ?h°pprT“1  ”d  a  !hon 

electron  pulse  (beam  energy,  750kV;  current,  100A;  radius  =  2mm-  length  =  2fYl^^  I  pP^=th^ugh  amuch  longer 
the  example  of  pulse  propagation  in  a  constant  parameter  undulator  where  the  Gatesian  n°  This  B  ™  contrast  to 

identified  with  sidebanding.  Variation  of  initial  pulse  width  shows  convergence  b  a  h  ^  UP  ““  ^S^ties 

the  slippage  of  the  radiation  pulse  through  the  electron  milse  the  50psec wide  output  pulse.  Because  of 

three  times  the  kinetic  energy  density  of  the  eXmonK ’ thC  ^  micr0wave  Pdse  intensity, -SGW/cntf,  is  about 

Keywords:  short  pulse,  millimeter  waves,  free  electron  laser,  high  power 


This  paper  studies  a  numerical  model  of  a  traveling  wave  hi  ah  M;„  __T 

optimal  efficiency  using  a  variable-parameter  (“tapered”)  undulator’and  wK’tSh^  EEL  whicnh  operaies  at  nearly 
pulse  with  a  smooth  spectrum  almost  free  of  sidebands.  The  hardware  would  inchX  ^  “}f,ense  ”cIean"  output  spike 
suitable  pulse  having  a  Gaussian  shape,  as  input  to  a  high  efficiency  FFT  trover  seed  source  which  supplies  a 

tapered  undulator.  Ourfindings  are  that  one  nSght 

-o  GW/cm-  and  FWHM  -50  psec  using  a  0.75MeV,  100 A  elLtrotTbXm  S  ™  57  P  Se  havn«  ^  P°we‘ 

wavelengths  wide,  suggests  potential  application  in  the  area  of  impuS  ^  narrow  pulse>  about  ten 

Under  certain  conditions,  the  free  electron  laser  IEEL1  oscillator  hoc  r.  . 
chaonc  high  power  "spike"  pulses  of  radiation  characterized  bva^Se^-^  found  10  Pr°vide  an  output  of  narrow, 
operation  of  a  EEL  oscillator,  experiment  [4]  as  well  as  numericahhm  h-  . 8X1  ^  sP®ctrum  [1-3].  Furthermore,  in  the 
regime  [5],  shows  that  the  EEL  can  operate  m  a  mode  characteSld  **  f31^  wel1  in*°  1116  nonlinear 

m  a  mode  that  has  higher  efficiency  and  a  wide  spectrum  The  latter  has  to  do  3  n3IT0W  sPectrum,  or 

been  observed  experimentally  [7  8]  and  which  is  also  found  ;n  ™  ^  ^  11x6  Sldeband  instability  [6]  which  has 

both  [14]  arise  from  slippage.  However,  *ere is -13]  since 
appropriately  chosen  taper  of  the  undulator  f  15-171  We  therefore  •  1  instability  can  be  stabilized  with  an 

both  a  high  power  pulse  as  well  as  SeSlS ’ ^  US£  °f  3  ^  -d^or  to  £do£ 

....  We  now  devel°P  a  numerical  model  which  establishes  how  such  a  FFT  nnicp 
milhmeter  wave  pulse  which  is  propagating  along  a  much  longer  rmico  of  «i  LL  Ph  033  ^  P^P^ed;  we  study  a  short 
PEL  resonance,  as  she  light  one ^  tl8)'  Al 

wavelength  X, .  We  shall  study  the  case  where  the  electron  ZLS!;  I  F  of  the  el®otrons  by  one  optical 

=  Nw  Xs,  so  that  essentially  no  radiation  ap^aro ataSrf 7  ***? **  °veraU  sUpPage <*«« 4 
have  considered  the  wealth  of  ^  (A  nUmber  of  rece"t  publications 

mvoive  "supenadiance"  (9-13).)  We  begin  by  spying  a  set  oflD 


(x,  y) 
ax 


d6j  (x,  y)  _  _i 

dx  ~ 


2PT^A  ssinyj 

TjPjll 

1  _  ^s(  1  —  Pill) 
“P  L  kwPji!  J 


(1) 


(2) 


250/SPIE  Vo  I.  2557 


The  above  equations  are  derived  directly  from  the  orifrirwl  in  .  T ,  , 

transforming  the  variables  z  and  t  into  new  independent  var4le? TU?  y^te^L^  M  by 

(Z  "  VJl-t)/lc  ;  here>  !c  =  V4^p  is  the  cooperation  length  which  is  defined  as  the  minimum  distli’ce^ 

which  an  electron  may  interact  cooperatively  with  the  radiation  ITCH  o  =  —  I  .  ,  _. 

F  Tr  v4ckw  )  ’  12  the  Piercs  parameter^  J3j[  = 

[l  -  (P2  -  2awascos\jlj)/7j ] m  is  the  axial  velocity  of  the  jth  electron,  u2  =  1  +  £  +  aj  ,  aw  and  a$  are  the 
normalized  vector  potentials  of  the  undulator  (eB^w  me2  )  and  radiation  field  (eEs/ksmc2);  r  is  the  rela  ’  '  • 
factor  of  the  jth  electron  ,  \j/j  =  6j  +  (b  is  the  relative  phase  of  the  jth  electron  with  respect  to  the  radiation  pulse^O^y) 

=  Ase1®  is  the  complex  amplitude  of  dte  radiadon  poise  with  As  =  dp  andp  is  the  phase  shift  of  dte 

radiation  poise.  The  angular  brackets  in  the  right-hand  side  of  equation  (3)  indSe  an  ensemble  average  over  ail  electro 
Fcr  the  other  quantities:  \  is  the  resonant  energy  of  electron  in  units  of  me?,  kw  =  2tr/l„  is  the  wave  nnmber  of  2 
undulator,  lw  is  the  undulator  period;  ks  =  2rtAs  =  ajc  is  the  wave  number  of  the  tadiadon  ptdse.  ^  is  the  radiation 


wavelength,  lw  and  A$  satisfy  the  resonance  condition  = 


^  w 


0+0 


COp  — 


-  (^2^efa 


m 


is  the  plasma 


frequency  for  ne  e.ectrons/cm  .  The  ongmal  FEL  wave  equations  have  been  obtained  assuming  the  field  amnlimde  ,•« 
^owly  varying  funcaon  of  time;  this  approximation  has  been  verified  in  the  application  for  the  short  n„i~  P 
discussed  here.  The  operation  of  the  FEL  described  here  is  on  the  borderline  tenveen  the  Cnmm™  Jo  pr0paSaGon 
and  if  an  actual  device  is  to  be  constructed,  the  beam  space  charge  effects  should  be  included  to  Raman  reSimesr 

for  this  illustrative  example,  we  have  not  incorporated  this  complicatS  nor  t£at  ^  ov^^T  ^  h°WeVer’ 

element  to  contain  the  waves;  both  items  are  comparatively  straightforward  to  include  in  the  ahn  V®£U1  e<  an  essential 
—op  wc  shad  own  ft.  ouipui  wave  JL,. 


oc  M  J!1®  numerical  simulation  is  based  on  the  computational  model  described  by  equations  (1)  to  (31  The  FFT  wnrW 
as  a  traveling  wave  amplifier.  The  spatial  distribution  of  simulation  electrons  has  a  rectaneular 
akcu  to be monoenergedc,  and  at  the  undulator  entrance,  they  are  uniformly  distributed  insidfdre  Ln  Mu  wim„T 
thousand  simulation  electrons  per  radiation  wavelength.  For  each  wavelength-size  "strip"  of  electrons  thfmi+wt, 
location  of  the  electrons  with  respect  to  the  radiation  field  is  uniformly  distributed  from  -n  to  K  The  Lnt^lZ 
is  "seeded"  inside  the  electron  beam,  with  its  initial  amplitude  a.cn  =  ICO  (9  kW/cm2  at  1  5mm  «»  J  ^  puise 
prose  and  width  can  be  varied  so  mat  we  can  study  the  evolution of  7„t2™ 

computational  results  the  electron  beam  puise  and  the  radiadon  spike  are  ploded  as  the  funedon  of  indeSnien  x 

and  y  respeedvely,  which  implies  two  moving  -windows’  with  the  former  at  the  speed  of  light  c  and  thflaSTat 
of  electrons  v„.  both  are  scaled  in  the  units  of  the  radiadon  wavelength  ks  from  their  leading  edge  Stu  to 
quanddes  are  tecoroed  at  various  undulator  positions,  these  results  will  describe  the  dme  evolution  of  Routes  m-evL 
mst  ntns,  we  have  earned  out  sunuladons  for  different  FEL  paramemrs,  including  different  lengths  of  to  S 
radmuon  puls«.  Our  results  for  -soperradianf  pulse  evolndon  are  in  a  ve*  good  agreement  witf,  dtTof 
authors[9]  and  provide  a  calibrauon-of  our  code.  at  or  previ0lls 

_  ,  _Pur  maif  ^erest  1S  to  investigate  an  initial  shon  radiation  pulse  propagating  through  a  sufficiently  lone  electron 
pu  se.  The  initial  radiation  pulse  has  a  Gaussian  profile  and  is  injected  into  the  rear  pan  of  the  electron  beam  pulse;  the 


peak  of  the  initial  pulse  is  located  at  y  =  170  Xs.  The  radiation  pulse  starts  from  that  position  at  the  undulator  entrance,  and 
then  moves  toward  the  front  of  the  electron  pulse  as  it  moves  along  the  undulator.  The  Fourier-transformed  spectrum  is  also 
Gaussian  and  has  a  central  frequency  ois=U6x  1012  sec-1.  Other  simulation  parameters  of  the  EEL  amplifier  are  listed 
in  Table  I,  where  a  representative  electron  beam  pulse  length  of  30cm  (-lnsec)  is  taken;  this  beam  pulse  is  200 
wavelengths  long,  and  so  the  optical  pulse  moves  only  halfway  through  as  it  traverses  100  undulator  periods. 


Table  I.  Simulation  parameters  of  radiation  spike 
propagating  in  electron  beam  pulse 


Electron  beam  energy 
Electron  beam  current 
Beam  intensity 
Electron  beam  radius 
Beam  pulse  length 


v=  2  5 
lb  =  100  (A) 
0.63GW/cm2 
=  0.2  (cm) 

Lb  =  200  X$  (30cm) 


Undulator  period 
Undulator  taper 
Undulator  parameter 
Undulator  length 


lw  =  1.5  -  1.1  (cm);  linear  ramp 
r|  =  2  xlO  (cm"1) 
aw  =  0.2  (constant) 


Radiation  wavelength 
Spike  length(FWHM,  intensity) 

Initial  pulse  amplitude 
Peak  spike  amplitude 


Xg  =1.5  (mm) 

Lr  =10Xg 

asO  =  10"4  (9  kW/cm2) 
as  =  0.06  OGW/cm2) 


Pierce  parameter 
Cooperation  length 


p  =  0.023 
lc  =  4  Xg 


Figure  1  shows  the  pulse  evolution  in  a 
constant  period  undularor  (lw  =  1.5cm.  at  four 
undulator  positions  Nw  =  0,  25,  50  and  100. 
In  these  figures,  the  electron  envelope  is  at  rest 
and  the  optical  spike  propagates  from  the  right 
to  left  as  Nw  increases.  The  radiation  spike 
retains  the  initial  profile  until  approximately  50 
undulator  periods  where  the  growth  becomes 
saturated.  Saturation  occurs  when  the  loss  of 
beam  energy  causes  the  interaction  to  drop  out 
of  resonance.  Examining  the  spectrum  of  the 
pulse  at  Nw  =  100  (see  right  side  of  Figure  3), 
we  find  the  spectrum  has  developed  prominant 
sidebands;  these  become  comparable  to  the 
carrier  power  by  150  periods. 

If  the  undulator  is  tapered,  we  expect  the 
growth  of  the  sideband  can  be  suppressed,  and 
in  Figure  2  we  show  the  results.  The  undulator 
period  is  linearly  tapered  froml.5cm  to  1.05cm 
in  150  periods.  All  other  conditions  are  the 
same  as  in  Figure  1  and  are  obtained  from 
Table  I.  In  this  numerical  experiment,  we  did 
not  optimize  the  undulator  taper  to  pursue  the 
highest  efficiency  enhancement:  instead,  the 
appropriate  taper  was  chosen  for  the  purposes 
of  obtaining  the  "cleanest"  pulse  (however,  this 
optimized  taper  is  not  very  different  from  the 
taper  which  extracts  maximum  energy  from 
the  beam) .  Unlike  the  untapered  case,  the 
radiation  pulse  at  high  amplitude  does  not 


radically  change  its  shape,  but  displays  a 
narrow,  regular  profile  throughout  the  slippage  region.  The  width  of  the  initial  Gaussian  pulse  is  FWHM  ~10XS  in 
intensity;  it  grows  in  amplitude  and  broadens  during  the  first  25  periods;  it  continues  to  grow  while  keeping  a  nearly  regular 
profile  as  it  narrows  in  width  between  50  and  100  periods.  We  have  tried  various  widths  of  the  initial  Gaussian  pulse,  as 
well  as  a  different  initial  amplitude  profile  varying  as  ~  l/cosh(ay)  (a  solitary  wave  solution  [19]),  and  we  find  the 
variation  of  initial  choices  converges  to  a  similar  output  pulse  profile  and  width,  although  the  number  of  undulator  periods 
needed  to  form  the  self-similar  pulse  profile  may  vary.  This  suggests  that  a  short  optical  pulse  may  eventually  evolve  to  a 
self-similar  spike  in  the  slippage  region  of  the  tapered  undulator.  The  characteristic  width  of  the  spike  is  similar  to  the 
prediction  of  the  Ginzburg-Landau  solitary  wave  theory  [20].  In  the  constant  period  undulator,  the  self-similar  feature  is 
lost  when  the  FEL  system  goes  out  of  resonance.  Returning  to  Figure  3  Oeft  side),  we  find  not  only  a  cleaner  spectrum,  but 


evidence  of  sideband  suppression  as  well. 


In  Figure  4  we  compare  the  cases  of  an  untapered  undulator  (dotted  line),  a  taper  as  given  in  the  Table  (solid  line), 
and  a  steeper  taper  (dot-dash  line),  about  20%  larger  than  that  given  in  the  Table.  Examining  the  electron  energy  profile, 
we  find  considerable  efficiency  enhancement  in  the  tapered  undulators;  however,  the  optimum  taper  yields  not  only  a  better 
pulse  shape,  but  a  higher  peak  power.  The  maximum  value  of  the  normalized  field  amplitude  a$  =  0.06  (~3GW/cm~) 
corresponds  to  an  intensity  higher  than  the  beam  kinetic  energy  intensity;  this  enhancement  in  power  is  caused  by  the 
slippage  of  the  radiation  pulse  over  “new”  electrons  as  it  moves  down  the  undulator.  The  ragged  profile  of  the  electron 
energy  in  the  constant  undulator  is  found  to  smooth  out  in  the  tapered  undulator.  The  strong  intensity  of  the  spike  forms  a 
very  deep  potential  well  which  may  trap  most  of  the  electrons  even  though  they  have  an  energy  spread. 
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The  overall  power  gain  of  the  spike  pulse  is -50dB.  This  raises  the  auestion  of  nr.jM  •  , 

noise,  and  how  it  affects  the  amplified  pulse.  We  have  introduced  particle  noise  1221  inro  rh  part3cularly  electron  random 
cause,  under  certain  conditions,  a  “spiky”  noise  signal  that  both  leads  and  trails  the  main  n»i<V  o' 3d{!1and  ^  that  this  can 
wings  of  the  Gaussian  pulse  can  be  amplified  under  noisy  conditions  to  high  amplitude  and  if  enerf?  ^sent  in  the 
noisy  energy  that  precedes  the  main  pulse  will  broaden  the  energy  distribution  of  the  electrons  that  rhiViff506-1?,1  gain* the 
resulting  in  degraded  amplifier  performance.  However,  starting  from  our  choice  of  initial  power ^  int0’ 
panicle  noise  has  no  effect  on  the  output;  but  this  is  not  true  at  much  lower  input  levels  and  the.  we  fuid  11131 

important  at  shorter  wavelengths  (we  have  found  that  a  steeper  taper  is  useful  to  inhibit  noise  growth  orVh°  becomes  more 
The  type  of  noise  we  have  introduced  into  the  simulation  is  Gaussian,  and  it  should  be  nnJ^Vrhsr  ,°n  wavelength). 
practice  may  have  particle  noise  in  excess  of  this.  *  noted  0121  3X1  electr°n  beam  in 

In  conclusion,  we  have  found  that  injection  of  a  single  short  radiation  pulse  into  a  Inner  , 

mput  of  a  tapered  undulator  traveling  wave  EEL  amplifier  should  result  in  the* development 

comparatively  "clean"  in  both  time  and  frequency  domains,  and  which  propagates  in  a  nearlv  self  sim'iar  ^  6  is 
undulator.  Because  of  the  slippage,  the  peak  pulse  intensity  is  enhanSdTay  excS  L“l^J^ 

Smce  the  pulse  spectrum  is  regular,  the  output  pulse  from  the  FEL  should  be  useful  for  techfocV appSons 
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ABSTRACT 

One  or  more  electron  bunches  passing  along  the  axis  of  a  dielectric-lined  cylindrical 
waveguide  are  shown  to  emit  picosecond  pulses  of  high  power  broadband  microwave 
radiation.  The  bunches  can  be  generated  by  an  S-band  rf  gun,  and  thus  spaced  from  one 
another  by  10.5  cm  in  a  macropulse  sequence;  or  a  single  more  intense  bunch  can  be 
generated  using  a  laser-ffluminated  photocathode  in  the  rf  gun.  Theory  is  developed  for  the 
excitation  of  modes  of  this  waveguide  which  propagate  a.  the  bunch  velocity,  from 
Cerenkov  radiation  emitted  by  a  single  intense  bunch.  A  train  of  psec  coherent  wake-field 
pulses  are  shown  to  follow  the  bunch,  when  the  waveguide  modes  have  nearly  constant 
spacing  in  frequency.  An  example  is  shown  for  an  alumina-lined  waveguide,  with  10  nC, 
3-15  psec  bunches  having  an  initial  energy  of  6  MeV.  Computadons  are  presented  of  the 
mode  spectrum  of  the  radiation  and  its  time  structure.  It  is  also  shown  that  measurements 
of  the  mode  spectrum,  or  of  the  energy  loss  of  the  bunch,  can  be  used  to  infer  the  axial 

density  profile  of  the  bunch.  Certain  features  of  the  theoiy  are  compared  with  the  results  of 
a  preliminary  experiment. 


1 


I.  Introduction 


If  high-energy  electrons  pass  along  a  channel  traversing  a  dielectric  medium, 
coherent  millimeter  and  microwave  radiation  can  be  produced  by  the  Cerenkov  effect,  if  the 
electron  speed  is  greater  than  the  phase  velocity  of  the  waves  [1, 2,3,4].  When  the 
electrons  are  in  a  continuous  beam,  linearized  analysis  for  the  fields  and  the  self-consistent 
motion  of  the  particles  leads  to  a  dispersion  relation  which  predicts  exponential  growth  or 
decay  of  the  waves.  For  the  nonlinear  regime  of  particle  trapping,  the  phase  relationship 
between  the  waves  and  the  particles  will  either  allow  the  extraction  of  energy  from  the 
electrons  and  emission  of  radiation,  or  the  acceleration  of  the  panicles  [5],  In  previous 
analyses  of  this  problem,  the  electron  stream  has  been  initially  unbunched,  and  the  radiation 
occurred  at  a  discrete  frequency  given  approximately  by  the  intersection  of  the  electron 
beam  dispersion  line  (co  =  (3ck?)  with  the  dielectric-loaded  waveguide  line.  However,  as 
we  shall  show  in  this  paper,  if  the  beam  consists  of  a  single  short  bunch  of  electrons,  or  a 
succession  of  such  bunches,  then  many  modes  of  the  structure  can  be  excited 
simultaneously.  By  proper  choice  of  waveguide  dimensions  and  dielectric  constant,  we 
find  that  these  modes  can  be  phased  so  as  to  constructively  interfere,  to  produce  a 
distinctive  sharply-pulsed  Cerenkov  wake  field  which  trails  the  particles  as  they  move 
along  the  axis  of  the  channel  in  the  dielectric  [6].  This  Cerenkov  wake  field,  which  here  is 
embodied  in  the  TM0m  modes  of  the  structure,  has  a  short-pulse  (psec)  time  signature  and 
MW-level  peak  power  for  available  experimental  conditions.  The  efficiency  of  this  high 
power  microwave  source  depends  on  the  waveguide  geometry,  length,  bunch  energy  and 
bunch  size.  In  the  example  discussed  below,  20%  of  the  electron  bunch  energy  is 
converted  into  radiation  for  the  parameters  chosen;  this  percentage  increases  as  the  length  of 
the  interaction  increases.  Broadband  psec  microwave  and  millimeter  wave  radiation  may 
have  applications  in  remote  sensing  and  spectroscopy.  Conversely,  measurements  of  the 
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radiated  spectrum  and/or  bunch  energy  loss  in  a  selected  dielectric  waveguide  segment  may 
be  used  to  infer  the  axial  bunch  length. 

In  what  follows,  we  analyze  a  system  depicted  in  Fig.  1  which  consists  of  the 
following  elements.  A  thermonic  or  photoelectric  cathode  rf  electron  gun  is  taken  as  the 
source  of  the  electron  bunch  or  bunches.  Such  an  rf  gun  is  typically  powered  by  a 
conventional  microwave  source  (e.g.,  a  klystron)  which  provides  5-10  MW  of  L-  or  S- 
band  power.  This  gun  will,  at  the  crest  of  the  rf  field,  emit  pulses  of  electrons  about  3-9 
psec  (1-3  mm)  in  length  which  are  spaced  apart  by  the  rf  period  (10.5  cm  for  the  S-band 
case).  The  energy  of  the  electrons  can  be  as  high  as  6  MeV  for  a  2-1/2  cell  rf  gun,  as 
confirmed  in  experimental  tests  in  the  Yale  Beam  Physics  Laboratory  [7].  The  S-Band 
source  can  be  pulsed  at  up  to  several  hundred  Hz,  with  a  macropulse  length  of  several 
microseconds.  The  stream  of  electron  pulses  provided  by  this  gun  (with  suitable 
downstream  focusing  elements)  is  directed  down  the  axis  of  a  hollow  dielectric  cylinder, 
with  outer  radius  R  and  inner  radius  a  ,  coated  on  its  outside  with  a  conducting  layer.  The 
waveguide  modes  of  this  structure  which  readily  interact  with  the  on-axis  electrons  are  of 
the  TM0m  class.  For  simplicity  in  what  follow-s,  we  shall  assume  the  waves  to  move  only 
along  the  device  in  the  direction  of  the  electron  flow;  e.g.,  there  is  no  counter-traveling 
wave  from  reflections.  Spent  electrons  can  be  defected  and  collected  at  the  end  of  the 
system.  The  axial  phase  velocity  of  the  waves  nearly  equals  the  speed  of  light,  so  the 
radiation  presumably  can  be  coupled  into  free  space  using  a  gentle  uptaper  in  the  inner  and 
outer  waveguide  radii.  However,  discussion  of  this  tapered  waveguide  and  the  resulting 
radiation  pattern  in  free  space  is  beyond  the  scope  of  this  paper. 

In  Section  II  we  present  the  wake  field  theory  in  such  a  cylindrical  waveguide  lined 
with  a  thick  shell  of  alumina.  Alumina  has  excellent  vacuum  and  mechanical  properties,  a 
high  relative  dielectric  constant  (-9-10)  and  low  loss.  It  is  a  ceramic  which  can  be 
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fabricated  into  the  geometry  we  have  chosen,  coated  on  its  interior  to  avoid  charge  build¬ 
up,  and  metalized  on  its  exterior  to  provide  a  good  outer  conducting  wall.  In  Section  EH  we 
present  computational  results  under  typical  conditions  for  a  single.injected  electron  bunch 
having  a  charge  of  10  nC,  that  will  be  seen  to  excite  fields  with  up  to  3.3  MW  peak  power. 
Wake  field  radial  and  axial  electric  field  patterns,  temporal  power  histories,  and  mode 
amplitudes  are  presented.  It  is  pointed  out  that  measurement  of  the  mode  spectrum  can  be 
used  to  infer  the  axial  bunch  profile,  and  that  measurement  of  the  energy  loss  of  a  bunch 
can  be  used  to  infer  the  effective  bunch  length.  In  Section  IV  we  describe  a  preliminary 
experiment  in  which  a  spectrum  of  modes  was  excited  in  a  dielectric-lined  cylindrical 
waveguide,  using  a  600  keV  unbunched  beam;  the  observed  mode  frequencies  are 
compared  with  frequencies  calculated  using  the  theory  given  in  Section  II.  In  Section  V 
results  of  the  work  are  summarized. 


II.  Wakefield  theory  and  radiation  modes 


The  model  analyzed  here  is  a  cylindrical  waveguide,  consisting  of  dielectric  material 
with  an  axisymmetric  hole;  the  outer  surface  of  the  cylinder  is  coated  with  a  low-loss 
conductor.  The  dielectric  constant  k  is  assumed  to  be  independent  of  frequency.  The 

geometry  is  depicted  in  Fig.  1.  Electrons  are  injected  along  the  z-axis  in  discrete  axi- 
symmetnc  bunches,  with  charge  density  given  by  p(r,z,t)  =  -(Ne/2nr)5(r)f(z-vt), 

where  e  is  the  magnitude  of  the  electron  charge;  N  is  the  total  charge  number  in  the  bunch; 
S(r)  is  the  transverse  charge  distribution,  assumed  to  be  of  infinitesimal  width  in  r;  and 
f(z  -  vt)  is  the  longitudinal  charge  distribution  for  bunch  particles  moving  at  axial  speed 
v.  For  this  geometry,  axisymmetric  orthonormal  wave  functions  can  be  found  for  the 
electromagnetic  fields  that  separate  into  TE  and  TM  classes.  Only  the  7M-modes  have  an 
axial  electric  field;  these  are  the  modes  considered  here.  Non-axisymmetric  modes  do  not 
separate  into  TE  and  TM  classes.  For  the  geometry  shown  in  Fig.  1,  conditions 
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can  be  found  where  all  TM  modes  have  phase  velocities  equal  to  v,  corresponding  to 
wakefields  that  move  in  synchronism  with  the  electron  bunches.  The  field  components  for 
the  complete  orthonormal  TM  mode  set  are  given  by 


where 


Ez(r,zj)  = 


m= 0 


/m(r)  -ia>n 

a„, 


?/v 


/*(0  = 


1  | Po(hm*R*a)  h{k\mr\  0<r<a 

P0{k2m,R,a)  [P0(k2m,R,r)  I0(klma),  a<r<R' 


(1) 


(2) 


and  z0  =  z-vr; 


Er(r,z,t)  =  i^jE„ 


m- 0 


g„,(r)  e-ia>mz0/v  ' 

a,n 


(3) 


where 


8m(r)  = 


1  f  7  P0(kZm,R,a)  Ii{kh„r),  0 <r<a 


P0(k2m,R,a )  [ /* .P, /0(/rlwfl),  a<r<R' 


(4) 


and 


Hg(r,z,t)  =  cy5  Er(r,z,t)- 


T,  0  <  r<  a 
K,  a  <r<  R 


(5) 


As  we  will  see  later,  the  field  amplitude  Em  is  expressed  by  the  product  of  the  Coulomb 
field  and  a  structure  factor  that  depends  on  the  electron  bunch  size.  In  the  above  equations, 
P0(k,R,r)  =  J0{kR)N0{kr)~  J0{kr)N0{kR)  and  P^R^)  =  J0(kR)Nl(kr)  -  J^N^kR)-, 
Jm(x )  and  Nm(x )  are  m-th  order  Bessel  functions  of  the  first  and  second  kinds,  and  Im(x) 
is  the  modified  Bessel  function;  a  and  R  are  radii  of  the  central  vacuum  hole  and  the  outer 
waveguide  wall,  respectively.  The  normalizing  constant  is  found  to  be 
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•“1/2  i/ 7 

where  /3  =  v/c;  y  =  (l-/32j  ;  and  =  (tc/3'  —  l)  ’.  The  (evanescent)  transverse 

wave  number  in  the  vacuum  is  kXm;  the  (real)  transverse  wave  number  in  the  dielectric  is 
k,m,  and  klm=C0m/cpy  =  k2mrK/r .  The  eigenfrequencies  are  a>m=cpyklm  =  cPrKkim. 
Since  all  TM  modes  have  phase  velocities  equal  to  v ,  we  also  have  com  =  c{3k.  m.  For  the 

fields  given  by  Eqs.  1-5,  ortho-normalization  obtains  in  the  form 


A 

I 


^rrE-m  ^z.n  & mn  f~  ~~~ 

y  am  an 


&0  k-n  *-xp[  *0  {.ZOnl  COn  )/v]. 


(6) 


where  D*zn  =  eE*.  n  -  K£0Ezn  in  the  dielectric,  and  D*.  n  =  e0E*zn  in  the  vacuum  hole. 
The  dispersion  relation  is  found  to  be 


Ix{kxa)  =  Kkx  JojkjR)^  {k2a)  -  7, feHM 
I0(kxa)  ~  k,  J0(k2R)NQ (k2aj-J0(k2a)JV0(k2R) 


(7) 


It  is  noted  that  one  can  have  eigenfrequencies  with  nearly  periodic  spacing,  since 
k2m(R  -  a)  ->  (n  +  l/2)zr  as  k  -»  °°.  As  the  asymptotic  eigenfrequency  spacing 

approaches  Aco  =  ncp^R- -1  .  The  wake  field  is  more  strongly  peaked  and 

more  closely  periodic  in  z0  as  the  eigenfrequencies  become  more  nearly  periodic,  i.e.  as  a 
higher  value  of  k  is  used. 


To  find  wake  fields  induced  by  an  electron  bunch,  one  expands  in  orthonormal 
modes  the  solution  of  the  inhomogeneous  wave  equation. 
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(8) 


A  <?2  Kir)  d2  1 

r  <9r 

l  <?rj 

+  dz1  c 2  dr  _ 

t)  =  S.(r,z,t)t 


withlhe  source  function  S.(r,z,r)  =  /i0%  +  — ' where  the  -component  of  the  current 

"  at  £a  oz 

density  is  jz(r,z,t )  =  vp(r,zj),  and  where  S;(r,z,r)  =  0  for  r  >  a.  A  complete  solution 
can  be  constructed  from  fields  as  given  in  Eq.  1-5,  since  these  are  solutions  of  Eq.  8  with 
S.(r,z,t)  =  0  everywhere.  We  expand  the  solution  of  Eq.  8  in  the  interval  0  <  r  <  R  in  a 

Fourier  integral. 

©o  00 

Ez(r,Z,t)=  £~fm(r)  jdkAm(l<)e~lkZo- 

_ A  M 


Inserting  Eq.  9  into  Eq.  8,  and  multiplying  both  sides  by  w(r)D.*„(r,z,r)  gives 


Am(k)  = 


2nam(k2  -(£>ljv2) 


jfdr'  jdz^Sz{iJrJ0)K(rVn,(r')»V)eikZ°  (10) 


where  the  weighting  factor  is  w(r)  =  [  -  y2,  fK]  in  the  intervals  [(0  <  r  <  a),  (, a<r<R )], 
respectively.  Then,  integrating  over  k,  with  due  regard  for  choice  of  the  contour  of 
integration  consistent  with  causality,  and  invoking  Eq.  6  yields 
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where  the  Green's  function  Gm(r 


(12) 


2  co,„ama 
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For  a  rectangular  bunch  p(r,z,t)  = -Ne(U  2nr)5(r)/Az  in  the  interval 
z0  -AzJ2  <  z<  z0  +&z/2,  and  p(r,z,t )  =  0  otherwise,  one  finds  for  the  spontaneous  wake 

field  the  result 


Ez(r'z'‘)  £*X  ’a 


f„,(r)  sin(m„,Az/2v)  io  . 


m  olv 


/»= 0 


(amAzJ2v) 


(13) 


While  for  a  Gaussian  bunch  p(r,z,r)  =  — - — 7- ^exp1 

2  TtrAz 


NeS{r)  \  (  ,  ^2 

- - —^-[-{z0/Az) 


,  one  finds 


I\  a,H 

n  1=0 


(14) 


In  Eqs.  13  and  14,  Ea  =  -Nej4n£0  a 2  is  the  Coulomb  field  of  the  bunch  at  the  edge  of  the 
hole,  and  causality  dictates  that  the  results  are  valid  only  behind  the  bunch,  i.e.  for  z0  ^  0; 
ahead  of  the  bunch  the  fields  are  of  course  zero. 


Eq.  7  has  been  evaluated  for  a  waveguide  with  a  =  0.128  cm,  R  =  1.0  cm,  and 
K—  10.0.  A  relative  dielectric  constant  of  k  =  10.0  is  close  to  the  value  of  9.6  for 
alumina.  It  is  assumed  that  K  is  independent  of  frequency.  The  phase  velocity  of  all  the 
TM  modes  corresponds  to  the  speed  of  an  electron  bunch  with  y  =  13,  i.e.  0.99704c. 
Fig.  2  shows  the  eigenmode  frequencies  for  this  choice  of  parameters.  For  this  case,  the 
first  frequency  interval  (cik  —coi)/2x  =  5.32  GHz,  while  the  asymptotic  interval  Aco/2x 

is  5.71  GHz.  It  is  thus  seen  that  conditions  can  be  found  where  the  eigenmodes  all  have 
equal  phase  velocities  with  nearly  equal  eigenfrequency  spacings.  The  maximum  deviation 
in  frequency  interval  for  the  example  shown  in  Fig.  2  is  6.8%,  about  double  the  value 
found  for  planar  geometry  [6]. 
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HI.  Power  output  and  field  structure 


The  theoretical  results  of  Section  II  were  evaluated  numerically  for  comparison  with 
wake-field  experiments  conducted  at  Argonne  National  Laboratory  [8].  Good  agreement 
was  obtained  for  the  wake  field  structure,  intensity  and  mode  frequencies.  In  this  case  only 
the  first  few  modes  make  significant  contributions,  as  also  noted  by  the  Argonne  group. 

The  detailed  computations  presented  in  this  section  are  mostly  for  a  Gaussian  bunch 
of  1.0  mm  length,  but  other  bunch  sizes  are  also  taken  for  evaluations  of  the  mode 
spectrum  and  energy  loss.  The  wake  field  structure  is  examined  for  a  total  charge  of  10 
nC,  although  the  results  can  be  scaled  in  direct  proportion  to  the  total  charge  (so  long  as  the 
bunch  energy  loss  is  not  so  great  as  to  significantly  reduce  the  bunch  velocity).  It  is  seen 
from  Fig.  2  that,  for  this  example,  consideration  of  only  the  first  few  modes  would  give  an 
incomplete  picture  of  the  wake  field  structure.  Therefore,  the  wakefield  is  computed  here 
by  including  modes  up  to  m  =  50  in  the  sum  in  Eq.  14,  even  though  modes  with  the  higher 
indices  will  make  smaller  contributions.  For  Q  =  10  nC,  the  Coulomb  field  E0  -  54.93 
MV/m.  The  computed  spontaneous  wake  field  pattern  of  a  single  bunch  is  shown  in  Fig.  3 
for  0  <  Zq  <  30  cm,  i.e.  at  the  instant  that  the  bunch  has  traveled  30  cm  into  the 

waveguide.  The  wake  field  peaks  are  seen  generally  to  alternate  in  sign;  to  each  be 
relatively  concentrated  in  z0 ;  and  to  have  a  period  of  10.5  cm,  corresponding  to  the 
vacuum  wavelength  at  2.856  GHz,  i.e.  half  the  asymptotic  frequency  interval  Aa/2n . 
The  peak  values  of  Ez  for  the  first  wake  at  the  bunch  and  at  5.2  cm  behind  the  bunch  are 
+20.51  and  -19.40  MV/m,  and  later  wakes  develop  secondary  oscillations  and  are 
somewhat  smaller.  The  waveguide  length  chosen  is  to  allow  a  moderate  20%  energy  loss 
from  the  bunch,  as  will  be  discussed  below.  The  radial  distribution  of  axial  electric  field  at 
different  axial  locations  within  the  first  half  wake  period  is  shown  in  Fig.  4.  This  figure 
shows  that  when  the  electron  bunch  passes  through  the  waveguide,  the  emitted  radiation 
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travels  in  a  pulse  moving  at  an  angle  with  the  waveguide  axis  until  it  reflects  from  the  wall 
and  heads  back  to  the  axis,  forming  a  radiation  cone  reminiscent  of  a  Cerenkov  cone  in  free 
space.  This  reflection  of  the  pulse  from  the  wall  back  to  the  axis  repeats  periodically,  and 
is  responsible  for  the  periodic  sharp  peaks  in  axial  wake  field  shown  in  Fig.  3. 


Fig.  5  shows  the  instantaneous  (time-dependent)  power  structure,  computed  as 


R 


P  =  2k  j  drr 


2 A 


V  n 


(15) 


It  is  noted  that  terms  with  m*n  contribute  to  the  instantaneous  power  and  are  responsible 
for  the  sharp  spikes  in  Fig.  5.  In  the  time-average  power  only  terms  with  m  =  n  contribute. 
For  this  example,  the  peak  power  is  3.3  MW  and  the  average  power  is  1.23  MW.  Fig.  5a 
shows  the  full  instantaneous  power  pulse,  which  occupies  a  time  duration  of  1  nsec.  Fig. 
5b  shows  an  expanded  plot  of  the  sharp  power  spike  near  the  center  of  the  pulse,  indicating 
the  3.3  MW  peak  with  a  width  of  3  psec.  Such  unusual  high-power  psec  pulses  could  find 
application  in  time-dependent  spectroscopy  or  remote  sensing. 

The  mode  spectrum  of  the  radiation  is  sensitive  to  the  bunch  length  and  profile. 
Shown  in  Fig.  6  is  the  wake  field  mode  spectrum  for  a  Gaussian  bunch  with  a  total  charge 
-Ne  =  Q  =  -10  nC,  for  several  bunch  lengths  Az  =  1,  2,  and  3  mm.  If  the  bunch  length  is 
increased  the  mode  spectrum  narrows,  with  diminished  amplitudes  for  the  higher  frequency 
modes.  Also,  the  narrower  the  bandwidth,  the  less  dielectric  dispersion  (neglected  here) 
will  matter.  The  strong  dependence  of  the  width  of  the  mode  spectrum  on  bunch  length 
could  be  the  basis  for  a  passive,  non-intercepting  diagnostic  tool  for  monitoring  bunch  size 
for  relativistic  electron  beams.  Profiles  with  non-Gaussian  shapes  will  have  different  mode 
spectra. 
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Now  we  turn  to  considerations  of  energy  loss  by  the  bunch.  The  rate  of 
energy  accumulation  with  distance  in  wake  fields  behind  the  bunch  dW/dz  is  given  by 
integrating  the  field  energy  density  over  the  waveguide  cross-section,  thus  leading  to  the 
following  relation. 


^  J d& J  rdr  |e(r) [e£/h  +  Ez.m  ]  +  VoH8.m  } 
m= 0  o  0 


e0EljcR2 

2 


I 

m=0 


h2{$m) 


at, 


(16) 


f 

+  K  ■ 
V 


Io{k\ma) 
A)(^2m*  R>a)  j 


R,r)rdr  i 


The  source  current  and  charge  distributions  in  E q.  8  are  assumed  to  remain  constant  during 
the  interaction.  If  not,  the  individual  mode  amplitudes  must  be  adjusted  iteratively;  this 
amounts  to  introduction  of  a  ^-dependent  structure  factor  h(%m)  in  equations  such  as  Eqs. 

13  and  14.  For  the  rectangular  bunch,  the  structure  factor  h(£m)  =  (sin<j;,rt/|m)2  ;  and  for 
the  Gaussian  bunch,  A(|/n)  =  ex p(-^,),  with  =  com As/2v .  For  the  parameters 

chosen  and  Q  =  10  nC,  Eq.  16  gives  dW/dz  =  4.01  x  10‘2  J/m.  Equating  dW/dz  to  the 
energy  loss  rate  QEdrag  gives  Edrag  =  4.01  MV/m  for  a  single  bunch.  This  drag  field 

will  slow  down  the  6  MeV  bunch  to  4.80  MeV  within  30  cm.  For  this  example  the 
efficiency  for  transfer  of  bunch  kinetic  energy  into  radiation  energy  is  seen  to  be  20%.  If 
the  waveguide  were  to  extend  to  100  cm  in  length,  this  drag  field  would  transfer  4  MeV  of 
beam  energy  into  wake  field  radiation,  giving  an  efficiency  of  67%.  Larger  energy  loss  by 
the  bunch  will  result  in  reduction  in  the  phase  velocities  of  the  wake  field  modes  that  are 
excited,  with  a  corresponding  change  in  the  wake  field  pattern. 
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Measurement  of  the  energy  loss  of  a  bunch  traversing  such  a  dielectric  waveguide 
can  also  be  used  as  a  bunch  length  diagnostic,  if  the  total  bunch  charge  is  known.  This  is 
illustrated  in  Fig.  7,  which  for  a  10  nC  Gaussian  bunch  shows  the  energy  loss  in  traversing 
a  30  cm  length  of  dielectric  waveguide  as  a  function  of  the  bunch  width  Ac.  It  is  seen  that 
energy  loss  is  a  strongly  decreasing  function  of  increasing  bunch  length.  This  diagnostic 
method  could  be  more  easily  implemented  than  the  measurement  of  radiation  spectrum, 
since  the  latter  requires  careful  calibration  of  the  detectors  over  a  wide  frequency  range, 
while  the  former  requires  knowledge  of  only  the  bunch  energy  before  and  after  traversing 
the  dielectric  waveguide.  Of  course,  the  spectrum  can  in  principal  be  inverted  to  cive  the 
actual  axial  bunch  profile,  rather  than  only  the  effective  width. 


IV.  An  experimental  test  of  certain  features  of  this  theory 

A  preliminary  experiment  was  carried  out  to  excite  Cerenkov  radiation  in  a 
cylindrical  tube  made  of  alumina,  using  an  unbunched  beam.  The  schematic  of  the 
experiment  is  shown  in  Fig.  8.  The  alumina  tube  has  outer  diameter  32  mm  and  fits  into  a 
copper  pipe  having  an  inner  diameter  of  33  mm;  the  axial  hole  through  the  alumina  cylinder 
is  9  mm  in  diameter.  The  alumina  cylinder  is  1 10  cm  long,  and  the  electron  beam  is  passed 
along  a  length  of  85  cm  before  it  is  deflected  to  the  wall  by  distorting  the  guide  magneuc 
field  using  an  iron  “beam  stop."  An  otherwise  uniform  guide  magnetic  field  of 
approximately  0.9  T  is  provided  for  electron  beam  guiding  and  focusing.  The  electron 
beam  has  a  diameter  of  about  3  mm  and  carries  a  current  of  about  100  A.  The  beam  is 
extracted  from  a  cold-cathode  diode  made  of  graphite,  which  fleld-emits  when  a  600  kV 
negative  pulse  is  applied  to  the  cathode.  Only  a  small  portion  of  the  emitted  current  is 
passed  through  an  aperture  in  the  anode,  a  procedure  which  has  been  shown  to  improve  the 
beam  quality  [9],  The  diode  is  attached  to  a  pulse-line  accelerator,  which  applies  a  150 
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nsec  pulse  of  high  voltage  to  the  cathode.  A  pressure  below  10  Torr  is  obtained  in  the 
diode  and  in  the  hole  through  which  the  beam  propagates.  Radiation  exits  the  waveguide 
into  air  through  a  plexiglas  window  which  is  cemented  onto  the  end  of  the  alumina 

cylinder. 

Radiation  was  initially  picked  up  by  several  detectors  connected  to  sections  of 
waveguide  used  as  high-pass  filters.  The  strongest  signals  were  obtained  at  S-band,  for 
which  the  TM01  mode  should  match  the  electron  speed  (0.97  c)  at  2.9  GHz.  From  the 
geometry  of  the  system  and  the  calibration  of  the  detectors,  it  is  estimated  that  MW  levels  of 
radiation  were  observed  for  pulses  lasting  <100  nsec.  However,  it  was  found  that  the 
radiated  power  could  be  propagated  in  X-band  and  K-band  waveguides  as  well.  This 
suggested  that  higher  modes  were  being  excited,  so  a  quasi-optical  grating  spectrometer 
was  installed  to  study  the  spectrum  of  the  radiation  at  wavelengths  shorter  than  12  mm 
[10].  An  example  of  this  spectrum  is  shown  in  Fig.  9.  The  horizontal  axis  is  labeled  in 
frequency  units,  with  the  frequencies  of  the  peaks  indicated;  the  spectrometer  was 
calibrated  against  a  24  GHz  source.  The  data  points  reflect  averages  with  a  standard 
deviation  of  roughly  30%.  In  addition  to  the  peak  frequencies  indicated  on  this  figure, 
another  peak  of  radiation  was  observed  at  39  GHz,  but  very  little  beyond  40  GHz.  The 
spacing  of  the  modes  is  seen  to  be  about  4  GHz. 

In  Table  I  we  compare  the  observed  peak  frequencies  radiated  versus  the  mode 
frequencies  predicted  from  the  dispersion  relation,  Eq.  7  of  Section  II,  using  a  dielectric 
coefficient  of  10  for  alumina.  The  frequency  spacing  predicted  was  about  4.3  GHz,  and 
the  mode  number  is  also  indicated.  The  theory  appears  to  be  in  relatively  good  agreement 
with  the  experimental  results  for  the  mode  frequencies. 


It  is  not  unreasonable  to  ask  why  this  system  should  radiate  high  power  at  the 
discrete  frequencies  as  predicted  by  theory,  since  the  beam  is  initially  unbunched  whereas 
the  theory  is  appropriate  to  either  one  bunch  or  to  a  sequence  of  bunches.  We  point  out  that 
the  observed  level  of  power  radiated  cannot  be  explained  without  the  presence  of 
considerable  bunching.  Evidently  a  linear  beam  Cerenkov  instability  is  present  that  leads  to 
rapid  bunching.  While  the  bunching  is  probably  not  nearly  so  distinct  as  used  in  the 
theory,  it  is.  apparently  enough  to  excite  a  mode  spectrum  up  to  m  -  9.  Thus  it  seems 
reasonable  to  compare  the  results  of  this  experiment  with  the  theory  of  Section  II;  but  we 
emphasize  that  this  experiment  does  not  check  other  features  of  this  theory,  such  as  the 
wakefield  or  the  time-structure  of  the  emitted  power. 


V.  Conclusions 


By  exciting  a  dielectric-loaded  waveguide  system  with  short  bunches  of  electrons, 
theory  predicts  that  high  power  microwave  pulses  of  Cerenkov  radiation  are  emitted  with 
picosecond-scale  time  variations.  This  occurs  because  the  short  bunches  excite  a  highly 
localized  wake  field  which  trails  the  electron  bunch  in  the  waveguide.  The  wake  field  is 
made  possible  by  the  coherent  superposition  of  many  discrete  TM  waveguide  modes.  We 
have  studied  the  example  of  6  MeV,  10  nC  bunch,  which  can  be  obtained  from  a 
conventional  rf  gun  driven  by  a  high  power  pulsed  klystron  at  2.856  GHz.  The  example 
showed  that  20%  of  the  kinetic  energy  in  the  pulse  can  be  converted  into  radiation,  with  up 
to  3.3  MW  peak  power  having  width  of  a  few  psec.  This  is  achieved  using  a  short  section 
of  waveguide,  and  only  a  single  bunch  is  analyzed.  When  more  than  one  bunch  is  present 
in  the  waveguide,  much  higher  radiation  power  may  be  reached.  To  our  knowledge,  this  is 
the  first  time  that  this  approach  to  the  generation  of  psec  MW-level  microwave  power  has 
been  discussed.  This  radiation  could  find  applications  in  material  studies  or  remote 
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Table  I.  Comparison  of  theoretical  and  experimental  mode  frequencies 


Mode  number 

ffi 

f  7 

h 

U 

Theoretical  (GHz) 

23.96 

28.32 

32.68 

37.06 

Experimental  (GHz) 

23 

27 

31 

35 
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FIGURE  CAPTIONS 


Fig.  1.  Schematic  for  a  Cerenkov  wakefield  wide  band  radiation  system.  The 

dielectric-lined  waveguide  is  used  for  wakefield  generation.  The  radii  of 
waveguide  and  empty  hole  are  R  and  a  respectively. 

Fig.  2.  Wakefield  mode  frequencies. 

Fig.  3.  Wakefield  radiation  from  a  10  nC,  6  MeV,  1  mm  long  Gaussian  bunch  after 

it  has  traveled  30  cm  left  to  right  in  the  cylindrical  waveguide. 

Fig.  4.  The  radial  distribution  of  axial  electric  field  at  several  axial  locations  of  the 

first  half  wake  period. 

Fig.  5.  (a)  Wakefield  power  distribution  from  a  10  nC,  6  MeV,  1  mm  long 

Gaussian  bunch  after  it  has  traveled  30  cm  left  to  right  in  the  cylindrical 
waveguide;  (b)  An  extended  plot  of  the  sharp  power  spike  near  the  center  of 
the  pulse. 

Fig.  6.  Wakefield  power  spectrum  for  variable  bunch  lengths. 

Fig.  7.  The  energy  loss  by  Cerenkov  radiation  as  a  function  of  bunch  length  for  a 

10  nC  Gaussian  bunch  traversing  a  30  cm  long  waveguide. 

Fig.  8.  The  schematic  of  the  experiment  demonstrating  the  wakefield  generation  in 

the  alumina-loaded  cylindrical  waveguide.  The  unbunched  beam  is  used  for 
this  experiment. 
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Fig.  9.  Power  spectrum  of  the  wakefield  emitted  from  the  experiment  shown  in 

Fig.  8. 
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Harmonic  millimeter  radiation  from  a  microwave  free-electron-Iaser  amplifier 
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Using  a  fires-eiectron  laser  iFEL)  configured  as  a  traveling-wave  amplifier,  we  have  caused  the  bunchin- 
produced  by  the  amplification  of  a  coherent  microwave  source  to  drive  appreciable  power  at  the  harmonics.  A 
lO-k'.V  24-GHz  microwave  input  signal  grows  to  the  -200  kW  level  using  the  tower  frequency  unstable  root 
of  the  waveguide  FEL  dispersion  relation.  The  FEL  operates  in  the  TEll  mode,  using  a  helical  undulator 
11.85-cm  period)  and  a  3-mm-diam  600-kV  electron  beam  contained  in  an  8.7-mm-i.d.  cylindrical  waveauide. 
Tne  harmonic  currents  set  up  by  the  microwave  afe  found  to  cause  growth  of  harmonic  power  under  two 
conditions.  First,  if  the  design  is  such  that  the  upper  frequency  root  corresponds  to  the  third  harmonic,  we  see 
small  amounts  of  third-  and  second-harmonic  power,  coherent  with  the  source.  Second,  we  have  found  kYV 
emission  of  the  seventh  harmonic,  most  likely  from  the  TE72  mode,  which  travels  at  the  same  speed  as  the 
2d-GHz  wave.  In  order  to  excite  the  seventh-harmonic  radiation,  the  electron  beam  must  be  displaced  from  the 
axis  of  the  guide  by  -  2  mm.  In  both  cases,  no  harmonic  power  is  produced  without  gain  at  the  fundamental. 
We  present  a  one-dimensional  theoretical  model  of  the  experiment,  and  use  the  numerical  results  to  interpret 
our  tindings.  The  mode:  predicts  that  if  the  microwave  signal  is  strong  enough  to  drive  the  FEL  into  saturation, 
the  harmonic  radiation  mould  become  powerful.  [S 1063-65  lX(97i  14108-3] 


PACS  numben si:  41.60. Cr 


INTRODUCTION 

If  an  electron  beam  undergoing  FEL  {free -electron  laser) 
interaction  is  bunched  by  a  coherent  source  of  power,  then 
FEL  harmonic  radiation  will  grow  from  an  initial  condition 
that  is  very  much  above  noise  input.  The  radiation  at  the 
harmonics  should  be  related  in  phase  and  frequency  to  the 
source,  which  bunches  the  electron  beam.  However,  lacking 
some  resonant  or  phase-matching  process,  typically  har¬ 
monic  radiation  remains  at  a  low  level  [l].  Appreciable  har¬ 
monic  power  depends  on  gain  at  the  FEL  harmonic  [2].  or  on 
seme  mechanism  to  couple  energy  from  che  lower  frequency 
FEL  wave  into  the  harmonic.  In  this  paper  we  study  an  ex¬ 
ample  of  each  of  these. 

Recently.  Pioveila  e:  at.  [3]  studied  a  waveguide  FEL 
where  the  beam  is  bunched  by  a  growing  wave  at  the  low 
frequency  beam-wave  intersection,  prompting  growth  at  a 
harmonically  related  high  frequency  upper  intersection  [see 
Eq.  11)  below]:  a  similar  idea  was  described  by  Stembach 
and  Ghaliia  [4].  In  the  waveguide  FEL.  there  are  two  un¬ 
stable  roots  of  the  dispersion  relation  that  represent,  in  the 
work  we  describe  in  this  paper,  growing  TEll  waves  at  24 
and  —72  GHz  (third  harmonic?:  the  frequencies  are  obtained 
approximately  by  calculating  the  intersections  between  the 
beam  wave  dispersion  and  the  elec  to  magnetic  wave  disper¬ 
sion  curves  in  a  waveguide,  neglecting  the  Raman  space 
charge  effect: 

(l) 

where  a»t-c0kw/(  1-/3)  and  kn.=l-/lw.  X  involves  the 
cutoff  frequency  of  the  waveguide.  X  =  [(l.84/ 
2—  R..)l„ //Jyjj'.  Here.  (3  c  =  v ;  the  axial  electron  velocity. 
'/;  is  the  corresponding  relativistic  factor.  lw  is  the  undulator 
period,  R3  is  the  cylindrical  waveguide  radius,  and  c  is  the 


speed  of  light.  Thus  one  might  generate  frequency  and  phase 
stable  millimeter  power  in  a  FEL  amplifier  using  a  micro- 
wave  source  for  bunching:  and  because  the  high  frequency 
harmonic  wave  in  the  FcL  is  also  grow*ins.  the  power  output 
in  principle  can  be  nigh.  V/e  have  reported  aireadv  on  an 
experiment  that  studies  such  harmonic  radiation  from  a 
FEL — configured  as  a  traveling-wave  amplifier — usine  a 
magnetron  source  to  drive  the  low  frequency  FEL  mode  [5]. 

In  this  paper,  we  give  further  data  and  interpretation  re¬ 
garding  the  experimental  properties  of  this  interaction.  .Also, 
our  recent  observation  [6]  of  appreciable  revenr/t-harmonic 
power  is  a  potentially  significant  development,  and  involves 
a  mechanism  that  is  different  [/]  from  that  described  above; 
we  shall  present  experimental  results  and  interpretation  of 
this  effect  as  well.  In  the  following,  we  present  a  theory  and 
a  numerical  simulation  of  the  experiments  that  makes  use  of 
the  FEL  equations  using  two  waves,  harmonically  related  in 
frequency  but  traveling  at  different  speeds  in  a  wavecuide. 

EXPERIMENT 

Figure  l  shows  a  schematic  of  the  FEL  apparatus,  which 
is  driven  by  a  pulseline  generator  [8],  This  pulseline  pro¬ 
duces  a  150-nsec.  nearly  fiat  voltage  pulse  (400-600  kVl. 


FIG.  1.  Schematic:  a  pulseline  source  provides  a  negative  150- 
nsec  pulse  to  a  graphite  diode  that  “apertures’’  the  electron  beam. 
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TABLE  I.  Parameters 

of  the  FEL. 

Col.  1 

Col.  2 

Electron  energy 

600  kV 

Low  frequency  FEL  mode 

24  GHz 

High  frequency  FEL  mode 

72  GHz 

Undulator  period 

2.5  cm 

1.85  cm 

aw,  effective  undulator  parameter 

0.2 

Guide  field 

0  G 

8800  G 

TEl  1  cutoff  frequency 

18.4  GHz  20.7  GHz 

Electron  beam  current 

-100  A 

Electron  beam  diam 

—  3  mm 

which  is  applied  to  a  cold  graphite  cathode  in  a  field- 
immersed  vacuum  diode.  The  solenoidal  field  is  primarily  for 
beam  transport,  but  it  also  slightly  enhances  the  spiralins 
motion  of  the  electrons  moving  through  the  undulator 
(quiver  velocity  —  0.1c).  The  electron  current  emitted  from 
the  cathode  is  “apertured"  by  a  3-mm-diam  hole  in  a  graph¬ 
ite  anode  to  improve  the  beam  quality:  the  beam  carries 
about  100  A.  Table  I  lists  parameters  of  the  device  appropri¬ 
ate  to  the  experiment,  which  is  conducted  in  the  ‘‘Raman** 
regime  (column  2).  whereas  the  first  column  refers  to  equiva¬ 
lent  parameters  for  the  ‘•Compton**  theory.  A  24-GHz 
10-kW  signal  obtained  from  a  magnetron  is  launched  in  the 
TE11  mode  inside  a  cylindrical  stainless  steel  wavecuide. 
enclosed  by  a  helical  unauiator.  A  relatively  small  value  of 
undulator  field  (  — 250G)  is  chosen  so  that  the  system  be¬ 
haves  like  a  traveling- wave  amplifier  with  a  gain  <50,  in 
which  case  oscillation  is  avoided:  also,  numerical  study  of 
the  electron  motion  has  shown  that  a  low  undulator  field  will 
give  stable  well-behaved  electron  orbits  in  the  combined  un¬ 
dulator  and  guide  fields.  The  propagation  and  alignment  of 
the  electron  beam  through  the  waveguide  were  studied  bv 
taking  a  series  of  “witness  plates/’  "(The  latter  determines 
the  beam  location  via  a  damage  pattern  on  a  piece  of 
thermal-sensitive  paper.)  The  diameter  of  the  waveguide  and 
undulator  period  are  chosen  so  that  there  is  exponential  2ain 
of  the  24-GHz  wrave  and  one  of  the  harmonics  at  an  electron 
beam  energy  of  600  kV:  a  single  frequency  code  [9]  that 
includes  the  one-dimensional  (ID)  motion  and  the  2D  wave 
equations  in  the  cylindrical  waveguide  was  used  for  this  de¬ 
sign.  Figure  2  shows  the  computed  gain  [9]  of  the  device  for 
beam  energies  of  600  and  400  kV,  taking  the  fundamental 
and  harmonic  waves  to  be  in  the  TEl  1  mode. 

Using  a  grating  spectrometer  adapted  for  millimeter-wave 
studies  [10],  we  have  found  24-GHz  power  gain  -20  (well 
below  saturation),  and  the  observation  of  the" first  three  har¬ 
monics  (—  100  \\ .  peak)  has  been  reported  [5].  No  harmonic 
power  is  obtained  unless  there  is  appreciable  gain  for  the 
24-GHz  wave  (compare  corresponding  traces.  Fig.  3).  The 
third-harmonic  power  output  was  found  to  be  greater  than 
the  second,  and  both  are  much  greater  than  the  fourth  (Fis.  4) 
or  higher  harmonics.  If  the  interaction  were  nonresonant, 
then  the  power  of  the  harmonics  should  fall  off  very  rapidly 
[1]  with  increasing  harmonic  number.  The  power  gain  at  72 
and  48  GHz  for  the  operating  condition  of  beam  voltaee  of 
600  kV  should  be  about  10.  according  to  Fig.  2;  this  figure 
shows  the  gain  maxima  that  result  from  the  two  intersections 
described  by  Eq.  (1).  Figure  3  shows  the  result  (lower  trace) 
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FIG.  2.  Computed  FEL  gain  vs  frequencv.  for  parameters  of 
Table  I.  column  2.  for  beam  energy  of  600  kV  (solid  line)  and  400 
kV  (dashed  line). 

of  mixing  the  72-GHz  FEL  third-harmonic  signal  with  the 
frequency  tripled  magnetron  signal,  demonstrating  a  fixed 
frequency  relationship  between  the  FEL  third-harmonic 
power  pulse  and  the  lower  frequency  magnetron  power  that 
drives  the  FEL  bunching. 

V^e  have  also  observed  the  emission  of  appreciable 
seventh-harmonic  power.  16S  GHz  (spectrum.  Fis.  4).  ob¬ 
tained  only  when  we  have  appreciable  gain  for  the  24-GHz 
FEL  wave  and  when  there  was  a  misalignment  of  our  system 
in  which  the  electron  beam  was  moved  about  2  mm  off  the 
waveguide  axis.  The  seventh-harmonic  signal  level  is 
- 1  kW  or  more,  whereas  the  24-GHz  signal  is  amplified  to 
200  kW.  much  less  than  saturation.  The  seventh  harmonic 
is  not  a  “resonant  FEL  interaction”  and  its  appearance 
should  depend  on  two  special  conditions.  The  first  is  that  the 


Time 

FIG.  3.  From  below:  bear  signal  envelope  of  FEL  third  har¬ 
monic  mixed  with  magnetron  tripled  frequency;  24-GHz  amplified 
signal:  diode  voltage  of  accelerator  (— 600kV);  72  GHz  detector 
signal.  Horizontal  scale:  50  nsec/div;  diode  voltage  vertical  scale: 
200  kV/div. 


56 


HARMONIC  MILLIMETER  RADIATION  FROM  A... 


2163 


Frequency  (GHz) 

FIG.  4.  Power  spectrum  scan,  covering  harmonics  2-8:  beam  is 
off  axis  for  harmonics  >5  (signals  for  the  seventh  harmonic  and 
harmonics  >5  are  zero  for  on-axis  beamf.  Data  points  represent 
average  of  several  shots,  recording  peak  power. 


temporal  and  spatial  spectrum  of  the  seventh-harmonic  elec¬ 
tromagnetic  wave  must  move  at  very  nearly  the  same  wave 
velocity  as  the  temporal  and  spatial  spectrum  of  its  electron 
current  source,  the  latter  being  set  up  by  the  seventh  tempo¬ 
ral  harmonic  of  the  bunching  caused  by  the  amplified  TE 1 1 
24-GHz  wave.  The  EM  wave  in  question  that  satisfies  this 
constraint  is  the  TE72  mode  of  the  cylindrical  waveguide 
[7],  since  this  wave  has  the  same  wave  refractive  index  as  the 
TE11  wave  at  24  GHz.  Also,  this  mode  has  zero  field  on 
axis,  and  therefore  the  electron  beam  must  be  moved  off  axis 
to  overlap  one  of  the  electric  held  maxima  of  the  TE72 
wave.  We  have  observed  the  seventh-harmonic  signal  in  this 
device  for  an  electron  beam  energy  of  600  kV  and  as  low  as 
400  kV  (Fig.  5),  since  there  is  gain  at  24  GHz  at  the  lower 
energy  (see  Fig.  2). 

In  contrast  to  the  *  ‘second-order* 4  FEL  resonant  wave  in¬ 
teraction  (which  commences  from  zero  source  current  and 
field,  apart  from  noise,  but  grows  exponentially),  the  har¬ 
monic  conversion  is  a  “first-order”  process  [11],  in  which 
the  beam  is  already  prepared  so  that  its  equilibrium  current 
carries  spatiotemporal  modulation  at  the  harmonic  (this 


modulation  is  set  up  by  the  FEL  resonant  interaction  at  24 
GHz).  Thus  there  is  no  need  for  a  “linearized  dispersion 
relation”  for  the  first-order  conversion.  What  must  happen, 
however,  is  that  the  wave  that  feeds  on  the  seventh-harmonic 
currents  in  the  beam  must  travel  at  the  same  speed  as  these 
currents  if  there  is  to  be  efficient  conversion.  The  FEL  inter¬ 
action,  which  generates  a  wide  spectrum  of  harmonics  as  the 
beam  becomes  progressively  more  bunched,  is  the  free  en¬ 
ergy  source  for  the  seventh  harmonic. 


THEORY 

To  study  the  coupling  between  the  harmonically  related 
waves,  we  solve  numerically  the  set  of  FEL  equations  with 
slippage  for  the  electron  dynamics  and  two  interacting  waves 
[5].  For.  simplicity,  we  disregard  the  space-charge  (Raman 
effect),  which  will  decrease  the  gain  and  require  a  shorter 
unduiator  period  to  provide  the  same  frequencies  as  used  in 
this  calculation.  The  parameters  of  the  calculation  are  given 
in  Table  I,  column  1.  The  low-frequency  (24  GHz)  TE11 
FEL  mode  has  group  velocity  —c/2,  as  does  the  seventh- 
harmonic  TE72  wave:  the  second-  and  third-harmonic  TE11 
FEL  wave  group  velocities  are  hisher  than  the  electron  speed 
(-0.88c). 

We  first  define  the  variables  .t  and  y  as 


x=A 


where  A  =  [ LH.( i/v « —  l/i?^)]*"1  and  v$x,vS2  are  the  group 
velocities  of  o)x  and  (x)Z*  is  the  electron  axial  speed,  and 
ZJ  are  the  axis  and  time  variables.  Denning  a{  and  a2  (ax 
=  as  the  normalized  vector  potentials  of  waves  1 

and  2.  which  have  axial  wave  numbers  kx  and  kz  in  the 
waveguide,  respectively,  and  aw  the  unduiator  potential 
eB.Jwl27rmcz ,  where  Bw  is  the  unduiator  field,  the  electron 
and  wave  equations  for  the  j th  electron  are 


lwk'tv2\  1 


(2) 


dv;  a,v 

77  =  -  U  [«!«., sin(a«:y+ d>y) 

T*  owi.T2$in(  0:;--r  <bZj)].  (3) 


da  i 
ay 


Time  ^ 

FIG.  5.  24-GHz  amplified  power  (below)  and  168-GHz  con¬ 
verted  power  (above),  obtained  when  the  diode  voltage  (middle)  is 
only  about  400  kV;  similar  result  can  be  obtained  at  600  kV  as  well. 
Same  scale  axes  as  Fig.  3. 


These  equations  also  involve  the  following  quandties: 
o>p  t  the  beam  plasma  frequency;  a.  the  harmonic  number, 
kwz  +  k2~  —  iuzt\  <j>X2i  the  opucal  phase;  and 

^2s^w+^"^2^t  Yri—k2il  +  a2w)I2kwl. 
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JIG.  6.  24-GHz  FEL  power  growth  (line  pair  above)  and  16S 
GHz  power  convened  (line  pair  below)  in  units  of  normalized  wav 
vector  potential  squared,  vs  axial  distance  measured  in  undulate 
penods.  There  is  a  1  %  mismatch  of  wave  group  velocity,  and  tw 
rases  are  shown:  for  initial  24-GHz  microwave  power  of  10  kV 
(broken  lines)  and  100  kW  (solid  lines). 

In  the  simulation,  the  FEL  behaves  as  a  traveling-wavi 
p  iner.  The  spatial  distribution  of  simulation  electrons  ha- 
a  rectangular  profile,  the  electrons  are  taken  to  be  monoem 
Seuc  and  at  the  undulator  entrance  ;  =  0  thev  are  uni- 
fonnly  distributed  inside  the  beam  length  with  100  simula¬ 
tion  electrons  per  undulator  period.  For  each  wavelength-size 
stnp  of  panicles,  the  relative  phase  location  of  the  electron* 
tth  respect  to  the  radiation  field  is  uniformly  distributed 
r-f?en  , 77  aund  “  For  *e  output  format  of  the  computa- 
*e  Jectron  pulse  and  the  radiation  pulse 
h,?aken  t0  be  fiat-topped)  are  plotted  as  a  function  of  the 
anables  x  and  y .  which  imply  two  moving  ‘  ‘windows.  ‘  ’  the 

rh^n^M  v*2  and. 1156  ,atter  at  ft  t  both  are  scaled  in  units  of 

f  axl  w?!0T  '  Th£  £leCtr0D  m°ti0n  is  FaralIel  t0  the 

x  axis  while  the  o>;  wave  moves  parallel  to  the  y  axis:  the 
u>)  wave  moves  along  the  characteristic 

XJ  1  ■  ./I  11" 

Wi  -Y=const‘ 

The  finite  radial  size  of  the  electron  beam  is  accounted  for  in 
the  code  using  a  multiplicative  “filling  factor”  for  the  radia- 
oon  current  term  on  the  right-hand  sides  of  Eqs.  (4)  and  (5). 

ile  this  can  be  estimated  from  geometry,  it  is  also  a  con¬ 
venient  parameter  to  adjust  the  gain  of  the  24-GHz  wave  so 
that  it  corresponds  to  our  observations. 

the^sevenrh  h  “  £Xampl£  *‘£  calculated  is  relevant  to 
tne  seventh-harmonic  experiment.  It  shows  the  peak  power 

grown  along  the  umiulaK)r  of  *e  24.  Md  168_g£ 

starting  from  an  initial  signal  of  10  kW  at  24  GHz  in  th* 
nght  circular  "rotating”  mode;  the  two  waves  travel  at 
near  y  the  same  speed  (the  group  velocity  is  about  0  5c  in  the 
waveguide).  The  figure  shows  L  the  conversion  of  r^ert 
wave  power  to  the  seventh  harmonic  is  -7%  when  thTri- 
reaches  the  end  of  our  undulator.  -32  periods  This  cal 
culation  (broken  line  pair)  is  for  a  wave^d  “ 


(a) 


FIG.  7.  Power  growth  of  24  (a),  and  72  GHz  (b)  waves  vs  axial 
distance  measured  in  undulator  periods.  The  three  examples  are  a> 
follows.  Solid  line:  initial  microwave  power  10  kW  and  train  of  40: 
dashed  line:  initial  microwave  power  100  kW  and  train  of  40:  dotted 
hne:  miual  microwave  power  10  kW  and  gain  of  400.  Tne  first  case 
is  that  of  the  experiment  where  - 100  W  at  72  GHz  is  produced:  in 
the  other  examples,  the  72-GHz  power  level  is  ~  10  kW  at  40  un¬ 
dulator  periods.  The  gain  of  the  24  GHz  wave  is  varied  bv  chancinc 
the  beam  filling  factor  0.2  for  low  gain  (experimental),  and  0.35  for 
the  higher  gain  example. 

(between  the  24-  and  the  168-GHz  signals)  of  1%.  approxi¬ 
mately  what  is  expected  from  the  theory  [7],  However,  if  the 
FEL  output  is  driven  near  saturation  (-2  MW.  solid  line 
pair),  the  computation  reveals  that  the  seventh-harmonic 
power  becomes  nearly  30%  of  the  microwave  power.  We 
point  out  that  this  is  a  ID  simulation  and  the  overlap  between 
the  particles  and  the  different  radiation  modes  is  being  mod¬ 
eled  with  the  same  filling  factor.  Since  a  3D  simulation  is 
really  necessary  to  model  the  off-axis  behavior  of  the  beam 
and  «s  overlap  with  the  waveguide  mode,  the  results  of  this 
calculation  cannot  be  compared  with  the  experiment  quanti- 


HARMONIC  MILLIMETER  RADIATION  FROM  A  .  . . 


2165 


rativeiy,  but  may  serve  as  an  incentive  for  further  study.  In¬ 
deed.  a  3D  code  [12]  has  been  used  recently  to  model  the 
generation  of  the  seventh-harmonic  signal.  Appreciable 
seventh- harmonic  output  (  — 2kW)  has  been  found  [13]  us¬ 
ing  the  Raman  version  of  this  code  with  the  electron  beam 
located  off  axis:  much  smaller  power  f  —200  W)  was  found 
when  the  beam  was  located  on  axis.  These  results,  which  are 
not  much  different  from  our  experimental  findings,  will  be 
presented  in  more  detail  elsewhere,  as  the  optimum  location 
of  the  electron  beam  for  power  conversion  has  yet  to  be 
found. 

Figure  7  shows  the  increase  of  24-GHz  ( power  gain 
=40)  and  72-GHz  power,  starting  from  10  and  100  kW  of 
the  24-GHz  signal.  The  two  waves  travel,  respectively,  at  a 
group  velocity  of  -c/2,  and  0.96c.  This  is  the  “beam-on- 
axis“  case,  where  the  panicles  couple  to  the  TEll  mode. 
The  third-harmonic  power  reaches  a  level  of  100  W  and  10 
kvV,  respectively.  Another  example  is  computed  where  the 
low  frequency  wave  starts 'off  at  10  kW\  but  experiences  a 
power  gain  —400.  By  comparing  these  examples,  one  finds 
that  in  order  to  get  harmonic  power  of  tens  of  kW  it  is  nec¬ 
essary  to  have  either  a  large  input  microwave  power  or  a 
high-gain  system:  the  requirement  is  that  the  device  achieves 
a  microwave  power  output  near  the  saturation  level  so  that 
the  harmonics  of  the  bunching  are  large  [1].  The  figure 
shows  that  the  harmonic  output  is  very  sensitive  tbvthe  output 
microwave  power:  we  find  that  as  much  as  70  kW  of  third 
harmonic  can  be  produced  if  the  computation  is  carried  fur¬ 
ther  along  to  50  undulator  periods.  The  simulation,  which 
compares  with  the  experiment  (solid  line),  shows  the  emis¬ 
sion  of  oniv  a  small  amount  of  power,  as  we  have  found.  The 
harmonic  power  begins  its  rapid  growth  when  the  bunching 
produced  by  the  microwaves  develops  a  substantial  harmonic 


component:  until  then  (at  Nw<20k  the  bunching  harmonic 
and  its  wave  amplitude  remain  small  and  “noisv1 2 3 * * * 7*  [[]. 

CONCLUSIONS 

We  believe  the  harmonic  power  observed  in  this  experi¬ 
ment  can  be  of  use  to  create  coherent,  phase-referenced 
power  in  the  millimeter  spectrum  for  applications  such  as 
radar  and  accelerator  physics.  Tne  higher  frequency  waves 
■  are  driven  by  the  bunching  set  up  by  the  microwave  source 
and  we  have  shown  in  the  example  of  the  third  harmonic  that 
this  signal  has  a  fixed  relationship  to  the  lower  frequencv 
source.  The  third-harmonic  radiation  does  not  compete  with 
the  seventh  for  the  free  energy  of  the  bunching,  because  the 
seventh  harmonic  does  not  occur  when  the  electron  beam  is 
positioned  on  the  axis.  On  the  other  hand,  one  can  suppress 
the  lower  harmonics  by  making  the  FEL  interaction  resonant 
at  one  microwave  frequency  only  [by  arranging  that  the  sec¬ 
ond  term  in  the  square  brackets  of  Eq.  i'I)  vanishes],  together 
with  positioning  the  electron  bean  q/f-axis  to  enhance  the 
seventh-harmonic  coupling.  Tne  latter  approach  may  repre¬ 
sent  a  very  convenient  way  to  generate  appreciable  2-mm 
wavelength  power,  using  a  FEL  operating  at  comparatively 
low  beam  energy  driven  by  a  powerful  microwave  source.  In 
our  experiment,  due  to  low  gain  of  the  system  (—13  dB)  and 
the  modest  microwave  power  available  at  the  FEL  input,  the 
power  output  of  the  harmonics  is  not  high.  However,  the 
theory  shows  that  one  can  expect  substantial  harmonic  out¬ 
put  if  the  FEL  amplifier  is  operated  so  that  the  input  micro- 
wave  signal  becomes  nearly  saturated  at  the  output. 
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